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In this paper we examine the accuracy of a variant of the depletion 
layer approximation as used in the analysis of the electrostatic potential, 
in the absence of mobile charge, in various charge coupled devices. The 
approximation is a linearization of the nonlinear potential equations 
and is important in the numerical solution of two-dimensional problems. 
A one-dimensional model is solved both exactly and by means of the depletion 
layer approximation and the two results are compared. We conclude that 
this variant of the depletion layer approximation is excellent for buried 
channel devices, and adequate though not as good for surface devices. 
Some criteria are given which indicate the excellence of the approximation 
and can be used to estimate errors in two-dimensional calculations. 

The purpose of this paper is to discuss the accuracy of a variant of 
the well known depletion layer approximation 1,2 used in the analysis 
of the electrostatic potential, in the absence of mobile charge, in various 
charge coupled devices (CCD's). 3,4 The approximation is a linearization 
of the nonlinear potential equations and is important in the numerical 
solution of two-dimensional problems. An application of these results 
is given in Refs. 5 and 6. We consider a one-dimensional model, which 
can be solved exactly, and examine in some detail the errors involved 
in the approximation. This analysis is a generalization of similar results 
in a previous paper. 1 

Consider a buried channel CCD 7 formed by a first layer, ^ y 5> h t , 
of Si0 2 ; a second layer, h t ^ y ^ h 2 , of completely ionized p-type 
Si; and a third layer, h 2 ^ y < °°, of uniformly doped Ji-type Si. (See 
Fig. 1) Then the dimensionless equations governing the potential 
<p(x, y) are 

VViCr, y) = 0, 

vW*i y) = <r(v)} 

VVsOc, y) = cxp (<p 2 (x, y)) - 1, 
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Fig. 1 — A schematic diagram of a buried channel CCD. 

where ip } (x, y) is the value of <p(x, y) in the j'-th layer and V 2 = d 2 /dx + 
d 2 /dy 2 . In eq. (1), the dimensionless potential <p(x, y), and distances, 
x, y, hi , h 2 , are related to the dimensional quantities (p*(x*, y*), x*, 
y*, h* u h* 2l respectively, (measured in MKS units) by 



<p(x,y) = -j-f<p*(x*,y*)> 



(2a) 



x = x*/\ D , y = y*/\ D , K = h*/X D , h 2 = h$/\ D , (2b) 
where the Debye length \ D is defined by 

X D - (eJcT/JNi)*, (2c) 

and e > is the electronic charge, k is Boltzmann's constant, T is 
the absolute temperature, e 2 the permittivity of the Si, and iV£ is the 
number density of the donors. Also, the dimensionless density of 
acceptor ions in the p-layer, <r(y), is related to the physical density 
Nt(y*) by 

a( y ) = Nt(y*)/NZ . (2d) 

We assume that at the boundary y = 0, <pi{x, 0) is completely specified. 
*fc0) = V (aO ^ 0. ( 3a ) 

while at the base y = °°, 
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tp>{x, oo) = 0. (3b) 

At the internal boundaries, y = fh , y = h 2 , <p satisfies standard electro- 
magnetic boundary conditions: 

ft(8, h,) = Vi {x, h), v ~ (•<-, /',) - ^ (*, /'.) + Q» , (4a) 

<p 2 (x, h 2 ) = <p 3 (.r, ft a ), XT (*i W = ^? (*i *0i ( 4b ) 

where 77 = 6i/e a is the ratio of the permittivities of Si0 2 to Si, and Q„ 
is a possible trapped charge at the oxide-semiconductor interface. The 
remaining boundary condition is taken care of by assuming periodicity 
in the x-direction 

<p(0, y) = tp(L, y), ^ y < cc (4c) 

(see Ref. 5 for a fuller discussion of these boundary conditions). 

In the buried channel CCD, the parameters are adjusted so that 
the potential minimum occurs in layer 2, h t 22 y ^ h 2 , where the mobile 
charges are to be trapped. For most purposes, it is only necessary to 
know the potential accurately in layer 2. The doping in the p-layer is 
typically introduced by ion implantation followed by drive in diffusion. 8 
Then the acceptor density can be given quite accurately by 

.(*,) = c. exp {-(^)* In c.}-l (5) 

where c, is the surface density of charge due to the original ion im- 
plantation. In some applications it is convenient to replace a(ij) by 
its average 

9 = y^-T I"' *W ( 'y = 2 -7T- erf ( V^) - 1. ( 6 ) 
"2 — >h J h t * vine, 

where erf is the error function. In a typical buried channel CCD, 
region 1, the oxide layer, is 0.1m thick and 77 = 1/3. Region 2, the 
p-layer is 5m thick with an initially implanted charge density at the 
surface of 4.6 X 10 I5 /cm 3 , which corresponds to an average doping 
of ~2 X 10 18 /cm 3 . Region 3, the n-layer, has a doping of 10 ,4 /cm 3 . 
In our dimensionless units, this corresponds to h x = 0.24, h 2 = 12.24, 
c, = 46 and a = 19.71 at room temperature. (In all the examples 
where a is used, we take a = 20 for simplicity). The potential has been 
scaled so that 1 volt corresponds to 40 dimensionless units. Typical 
plate voltages are in the range —400 % V„ £ 0. Although we shall 
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typically assume Q„ = in this paper, we include it because it is 
necessary in some applications. 

If layer 2 is removed, we have a surface CCD. 3 In a surface CCD, 
the mobile charge is trapped at the oxide-semiconductor interface, and 
only near this interface does the potential need to be known accurately. 
The remaining dimensions are typically the same as for a buried channel 
CCD, and the plate voltage is in the range -400 ^ V ^ -160. 

If cp 3 (x, h 2 ) is large and negative for ^ x ^ L, then for y > h 2 
and near h 2 , exp (<p 3 (x, y)) « 1 and the equation for <p 3 is approximately 
V V 3 = — 1. The region where this holds is called the depletion region, 
and at any point x, the depth of the boundary of this region below 
y = h 2 , R(x), is called the depletion depth (see Fig. 1). For (x, y) such 
that y > R{x) +h a ,\ <p 3 (x, y) | « 1 and we have approximately V V = 
<p 2 . In most problems R(x) is not a constant, but it is nevertheless 
often a good approximation to replace R (x) by some suitable constant 
R. This suggests replacing the system of eqs. (1) by the system of 
linear equations 

V 2 fa(x, V) = 0, < y < h, , (7a) 

S7 2 fa(x,y) = °(y), h<y<h 2 , (7b) 

V 2 + 3 (x, y) = -1, h 2 < y < h 3 = h 2 + R, (7c) 

V 2 +*(x, y) = Ux, y), h 3 = h 2 + R < y < <* , (7d) 

(See Fig. 2). In addition to ypi , fa and fa satisfying boundary conditions 
(3a), (4a), (4b) and (4c), we have the boundary conditions 

fa(x, h 3 ) = fa(x, h 3 ), ^ (x, h 3 ) = ^ (x, h 3 ), fa(x, oo) = (8) 

which must hold for ^ x ^ L. 

This paper is devoted to studying how the choice of the pseudo-depletion 
depth R affects the accuracy with which \p(x, y) approximates <p(x, y). 
We do this by considering the special case V (z) = V < 0, where V 
is constant, in boundary condition (3a). For this choice of V (x), eqs. 
(1) and (7) reduce to systems of ordinary differential equations the 
solutions of which are functions of y only. In all that follows, we drop 
all reference to x, replace V 2 by d 2 /dy 2 in eqs. (1) and (7), and consider 
only the one-dimensional problem. The remainder of the paper is 
devoted to calculating and comparing the values of ^2(^2) and ^2(^2) 
with ^2(^2) and ^2(^2), where ' = d/dy. It is easy to show that for 
h ^ y ^ h 2 , 
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Fig. 2 — A schematic diagram of the depletion layer approximation to a one- 
dimensional, buried channel CCD. 



\<p(y) - +(y) I ^ \v»(hj - 1>M I, 
<p'(y) - t'to) I ^ |rf(M - t'M I- 



(9) 



We only give the analysis for the buried channel CCD, that is very 
similar for the surface CCD. We give numerical results for both. 
We begin by solving for <p 2 (h 2 ) and <p' 2 {h 2 ). We can now write 

¥>i(y) = ay + V , (10) 

*(y) = r (S - v)<K0 ^ + &(i/ - h 2 ) + c. (11) 

J u 

In addition, a first integral of (lc) subject to (3b) is 1 

\WMY = exp (*(y)} - ^3(2/) - 1. (12) 

Making use of the boundary conditions (4) , it is now easy to show that 

a - - (b - JT*' crtt) d? + Q..) , (13) 

c = (jh - A, + ^)b + V + ^ Q.. - J[* - (* - fe, + ^)<r(0 ft, (14) 



ift 2 = e c - c - 1. 



(15) 
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Equations (14) and (15), for c = ip 2 (h 2 ) and b = <p' 2 (h 2 ), can be solved 
numerically. In fact, for most cases of interest c « — 1, so the term 
e c in (15) can be neglected, yielding 

b = <p'M = -(hi - h x +^) 

+ j(, 2 _ /ii+ y^ 2 [ Vo+ i-/;(^/ (1+ y^^ + ^.j , ( i6) 

c = v t (h t ) = -$b* - 1. (17) 

It can be shown that for a surface CCD, <^(/i,) and <p 2 (lh) can be ob- 
tained from (16) and (17) by setting o(y) = and h t = h 2 . 
To solve eqs. (7), we can write 

^(y) = Ay + V , ^ y ^ h, , (18a) 

Hv) = f' G ~ vMQ # + B(z/ - k) + c > >'' = y = h * ' (18b) 
My) - -Ky - k) 2 + Dfo - W + E > k ^ » ^ *■ - ( 18c ) 

M) = te h ' + *-\ h 3 = h 2 + R ^ i/ < * , (18d) 

where the boundary conditions at y = and y = «» have already 
been used. Upon employing the remaining boundary conditions at 
y = fh , h 2 and h 2 + R = fe a , we obtain six equations relating the seven 
constants A, B, C, D, E, 5 and R. 

Our main interest is in cases where R is given, so that A, B, C, D, E 
and 5 can be determined in terms of R. However, we first consider the 
case where 5 is specified. This will provide useful information about 
picking an optimal value of R later on, and will also give a useful 
criterion for estimating how good the approximate solution is when R 
is specified. 

In Ref. (1) the effects of requiring that 5 = (the usual depletion 
layer approximation) or 5 = - 1 were studied, and it was shown that 
the choice 5 = - 1 yields a better approximate solution. It should be 
noted that requiring 5 = —1 is equivalent to requiring that $"Qh) = 
WOk). We show here that if the value of 5 is specified, the choice 
5 = — 1 does yield the best approximation, but the values of \l/ 2 (h 2 ) 
and H(M are relatively insensitive to the choice of 5. It is easy to 
show that in terms of 5, 
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i; = • [i>2 - ih +-J -i 



lh - ih + ^ I + 5 J + 25 - 2V 

7}/ 



+ 2 f*' (t- h t + ^)<r($) dfcj - 2 ^ Q..]*, (19a) 
A = - R + - [Q.. - [' a® % - b) , 



^(/, 2 ) = B = D = R - 5, 

^(h,) = C = E=5 + R5- JR a . 



(19b) 

(19c) 
(19d) 



For a surface CCD, the expressions for yf/' 2 {lh) an d ^ 2 (/'i) are given 
by (19) when o(y) = and lh = lh • 

As a first example, pick the typical parameters mentioned earlier, 
/i, = 0.24, h 2 = 12.24, Q a , = 0, v = 1/3, V n = 0, and let <r(y) = d = 20. 
In Fig. 3, using eqs. (!(>)> (17) and (19), we show a plot of 



and 



«(*) = I Wh) - fP 2 (/« 2 ) |/| *■(&.) 



e(*') = | VM) - V ' a (h 2 ) |/| *> 2 (A 2 ) 



(20a) 



(20b) 



as functions of 5 for — 10 ^ 5 ^ 10. It can be shown from (17) that 
these relative errors have a minimum at 5 = — 1 (and are even functions 




Fig. 3 — Plots of c(\f/) and c(if>') for a buried channel device with lu = 0.24, 
hi = 12.24, ff = 20, Q ss = 0, 7? = \, and V = as functions of 5 for -10 ^ 5 ^ 10. 
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of 8 + 1), but over the whole range of 5 shown, the maximum relative 
error in \p 2 (h 2 ) is only 1.27 percent and in ^(/<„) is only 2.28 percent. 

As a second example, consider a buried channel CCD with the same 
parameters as in example one, except that in this case we assume that 
a{y) is given by eq. (5) with c. = 46. In Fig. 4, we show a plot of e($) 
for this case as a function of 5 for -10 ^ 8 ^ 10. The results are 
qualitatively the same as in the first example, except that the relative 
error is about twice as large. It should be remarked that the potential 
distributions in the first and second examples are quite different. The 
potential minimum in example one occurs much nearer the bottom of 
the p-layer than in example two. In example one, <p(h 2 ) = —1034.5 
while in example two, ip(h 2 ) = -530.6, so we would expect the depletion 
layer approximation to be better in example one than in example two, 
as it is. 

As a third example, consider a surface CCD with hi = h 2 = 0.24, 
a = 0, Q„ = 0, v = 1/3, and V = -160. In Figs. 5a and 5b, using 
eqs. (16), (17) and (19), we again show plots of e(f) and e(\p') as func- 
tions of 8 over — 10 ^ 8 ^ 10. We see a maximum relative error of 
1.53 percent in ^ 2 {h x ) and 18.3 percent in ^(^1) • 

Since the approximation is based on the assumption that <p(h 2 ) is 
large and negative, the magnitude of <p{h 2 ) must affect the accuracy 
of the approximation. For the buried channel CCD, <p{h 2 ) is < —500 
for all V ^ 0, so in this case, the approximation is very accurate for 




-12 -10 - 



Fig. 4 — Plot of c(\p) for a buried channel device with ft, = 0.24, hi = 12.24, 
c. = 46, Q„ = 0, v = h and V = as a function of S for -10 ^ 5 ^ 10. 
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Fig. 5(a) mid 5(b) — Plots of r(^) and <(\p'), respectively, for 11 surface device with 
h x = h 2 = 0.24, <r = 0, Q, s = (),, = • and V„ = -160 as functions of a for 
-10 < 5 g 10. 



any negative V . However, in the surface CCD, as V — > we have 
<p(lh) — » 0. Thus for "small" V , we expect the approximation to be 
quite bad. Figure 6 gives a plot of the relative error e(\p) for a surface 
CCD as a function of V using the same parameters as given for Fig. 5 
and 8 = — 10. To maintain an error of the order of 1 percent for all 
5 t [—10, 10], it is clearly necessary to keep V„ ^ —160. 

We turn now to the main problem where R is arbitrarily specified. 
This occurs typically in the numerical solution of two-dimensional 
problems where it can be reasonably assumed that the depletion depth 
R(x) does not vary too much with the various plate potentials. Then 
an average depletion depth is estimated, say from eq. (19a) by setting 
5 = — 1 and choosing for V„ some average plate potential. 

We now investigate the errors involved in this approximation. If 
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10-1 




-960 



Fig. 6 — A plot of c(\f>) for a surface device with ki = h% = 0.24, a = 0, Q, 
v = i and S = -10 as a function of V for -800 ^ V ^ -40. 

R is given, then it is easy to show that 

MM = b = d = [£' (*-*, + tya® « - ^ Q« - Vo 

+ iR 2 + ft]/[fi + 1 + fh - fh + ^] , 
A = - (b - P a(0 dH + Q, 



- o, 



fc0g = C = E = £R 2 + R - (R + 1)B, 
5 = R - B. 



(21a) 

(21b) 

(21c) 
(21d) 



As before, the above expressions give ^ 2 {hx) and ^(/?,) for a surface 
device when we set a = and /ii = h 2 . 

As a fourth example, we consider a buried channel CCD with hj = 
0.24, h 2 = 12.24, a(y) m a = 20, Q„ = 0, 77 = \, and we choose R = 
50.98. From (19a), this corresponds to the "true" depletion depth R 
(with d = -1) for a plate voltage V = -400. In Fig. 7, using eqs. 
(16), (17) and (21), we give plots of ef» and e(rp') as functions of V 
for — 800 ^ Vo ^ 0. In Fig. 8, we give a plot of 5 as a function of V 
from (21d) for the same data, and a plot of the "true" value of R 
from (19a) for 5 = - 1 as a function of V . Several important con- 
clusions can be drawn from these graphs. First, over -800 ^ V ^ 0, 
\l/(y) is an excellent approximation to <p(y) for hi ^ y ^ h 2 , the maxi- 
mum relative errors in \f/ 2 (h 2 ) and f' 2 (h 2 ), which occur at V = 0, being 
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Fig. 7 — Plots of e(\f>) and c{\p') for a buried channel device with hi = 0.24, 
h t = 12.24, a = 20, Q„ = 0, 77 = % and It = 50.98 as functions of V for -800 ^ 
V ^ 0. 

0.41 percent and 0.74 percent respectively. Second, ^{y) approximates 
ip{y) excellently over hi ^ y ^ h 2 , even though a relatively large error 
has been made in estimating the true value of R. In fact, R = 50.98 
differs from the "true" value by 11.6 percent at V = —800 and by 




d>(0)( = V ) 



Fig. 8 — Plots of 5 and "true" Ft for a buried channel device with hi = 0.24, 
hi = 12.24, a = 20, Q„ = 0, r, = J and l\ = 50.98 as functions of V„ for -800 ^ 
V ^ 0. 
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• (0) (=V ) 



-960 



Fig 9 — pi t, of c(^) for a buried channel device with h\ = 0.24, hi = 12.24, 
c. = 46, Q„ = 0, v = i and R = 39.66 as a function of V for -800 ^ V ^ 0. 

12.7 percent at V = 0. Finally, an estimate of the error in \f>(h 2 ) can 
be obtained from monitoring 5. For V ^ 0, as long as -7 ^ 8 ^ 5, 
the error in ^ 2 (/i 2 ) and M(h 2 ) is less than 1 percent. 

Next, consider a buried channel CCD with ft, = 0.24, h 2 = 12.24, 




<M0)(=V ) 



-960 



Fig. 10— Plots of 5 and "true" R for a buried channel device with /i x = 0.24, 
hi = 12.24, c. = 46, Q„ = 0, * - i and R = 39.66 as functions of V„ for -800 S 
V g 0. 
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Fig. 11(a) and 11(b) — Plots of ety) and c(\p'), respectively, for a surface device 
with hi = hi = 0.24, <r = 0, Q„ = 0, r, = $ and R = 26.55 as functions of V for 
-800 ^ V ^ -40. 



Q,, = 0, 77 = |, with <j{\j) given by (5) with c. = 46, and with R = 
39.66. From (19a), this corresponds to a "true" value of R (with 8 = — 1) 
for a plate voltage of V„ = —400. In Fig. 9, we plot e(^) as a function 
of V for —800 ^ V„ ^ 0. In Fig. 10, we give a plot of 8 as a function 
of V from (21d) for the same data, and a plot of the "true" value of R 
from (19a) for 8 = — 1 as a function of V . The results are qualitatively 
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the same as for example four, but as we should expect, the errors in 
using the depletion layer approximation are greater in example five 
than in example four. 

As a final example, we consider a surface CCD with /*, = h 2 = 0.24, 
a- = 0, Q„ = 0, ij = i and R = 26.55. From (19a), this corresponds 
to a "true" value of R (with 8 = - 1) for a plate voltage of V = -400. 
In Figs. 11a and lib, using eqs. (16), (17) and (21), we give plots of 
e(» and e(*') as functions of V for -800 ^ V ^ -40. In Fig. 12, 
we give a plot of 5 as a function of V„ from (2 Id) for the same data, 
and a plot of the "true" value of R from (19a) (for 8 = -1) as a func- 
tion of V . Hence, we note that over the range —800 <i V ^ —160, 
the error in ^ 2 (/'i) remains remarkably small, less than 1 percent. 
However, ^£(A,) does not approximate <p' 2 (Jh) nearly as well, especially 
as v — ► 0. Note also that in this case much larger errors have been 
made in estimating R, the error being 39 percent at V„ = - 160 and 
44.2 percent at V = -800. Again, 5 provides a check on the accuracy 
of the approximation, with - 10 < 8 < 8 assuring an error of less than 
1 percent in approximating <p(y) in /;, ^ y £ h 2 . 

We can conclude that this variant of depletion layer approximation 
works excellently for the buried channel devices in providing approxi- 
mations to <p(y) and if>'(y) in the region h t ^ y ^ h 2 . In addition, the 
parameter 8 is a good index of the excellence of the approximation. 
If in a two-dimensional calculation this variant of the depletion layer 
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*(0) ( = V ) 



Fig. 12 — Plots of S and "true" It for a surface device with h t = hi = 0.24, a = 0, 
Q„ = o, r, = J and R = 26.55 as functions of V for -800 ^ V ^ -40. 
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approximation is used, then a value of the pseudo-depletion depth R 
is estimated, say from using an average plate potential in eq. (19a). 
Then, so long as the solution of this approximate problem remains 
in the interval [— 10, 10] along the line y = R + lh , one can feel reason- 
ably confident that the approximate solution differs from the true 
solution, in the region of interest, by at most a few percent. It should 
be noted that the depletion layer approximation is poor near y = 
h 2 + R, but in many important cases this is unimportant. Finally, 
for surface CCDs, the approximation is good to tp{y) but not as good 
to <p'(y). 
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